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Abstract. We solve a family of Gaussian two-matrix models with rectangular 
N x (N + v) matrices, having real asymmetric matrix elements and depending on 
a non-Hermiticity parameter fi. Our model can be thought of as the chiral extension 
of the real Ginibrc ensemble, relevant for Dirac operators in the same symmetry class. 
It has the property that its eigenvalues are either real, purely imaginary, or come in 
complex conjugate eigenvalue pairs. The eigenvalue joint probability distribution for 
our model is explicitly computed, leading to a non-Gaussian distribution including K- 
Bessel functions. All n-point density correlation functions are expressed for finite N in 
terms of a Pfaffian form. This contains a kernel involving Laguerre polynomials in the 
complex plane as a building block which was previously computed by the authors. This 
kernel can be expressed in terms of the kernel for complex non-Hermitian matrices, 
generalising the known relation among ensembles of Hermitian random matrices. 
Compact expressions are given for the density at finite N as an example, as well as its 
microscopic large- N limits at the origin for fixed v at strong and weak non-Hermiticity. 



1. Introduction 

Non-Hermitian Random Matrix Theory (RMT) introduced by Ginibre [1] is almost 
as old as its Hermitian counterpart. At first it was seen as an academic exercise to 
drop the Hermiticity constraint and thus to allow for complex eigenvalues. However, in 
the past two decades we have seen many applications of such RMTs featuring complex 
eigenvalues precisely for physical reasons, and we refer to [2] for examples and references. 
Because matrices with real data are often modelled by RMT one could view the real 
Ginibre ensemble of asymmetric matrices as being the most interesting non-Hermitian 
ensemble. Unfortunately it has also turned out to be the most difficult one, as it 
took over 25 years to compute the joint distribution of its eigenvalues [3j 0], being 
real or coming in complex conjugate pairs. The integrable structure and all eigenvalue 
correlation functions were computed only very recently for the real Ginibre ensemble 

la ei ei ei ei ma E]. 

Our motivation for generalising this model is as follows. In the 1990's Verbaarschot 
proposed extending the three classical (and Hermitian) ensembles of Wigner and Dyson 
to so-called chiral RMT [12], in order to describe the low energy sector of Quantum 
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Chromodynamics (QCD) and related field theories. These chiral ensembles are also 
known as Wishart or Laguerre ensembles. Their non-Hermitian extensions [T3J, [T3] 
were motivated by adding a chemical potential for the quarks, which breaks the anti- 
Hermiticity of the Dirac operator in field theory. It was observed numerically quite 
early [H] that these chiral versions of the Ginibre ensembles have distinct features, 
either attracting eigenvalues to the real and imaginary axes (real matrices), repelling 
them (quaternion real matrices) or having no such symmetry (complex matrices). Only 
later was it realised how to solve these chiral non-Hermitian RMTs analytically, by using 
replicas [15] or by extending the initial one-matrix model plus a constant symmetry- 
breaking term [T3l [Ej to a two-matrix model. This idea from Osborn [16] led to 
a complex eigenvalue model that can be solved using orthogonal polynomials in the 
complex plane [17] . The solution of the two- matrix model was then derived for complex 
[T6l [T8] and quaternion real matrices [19J. Our paper aims to solve the third and most 
difficult of such non-Hermitian RMTs, a chiral two-matrix model of real asymmetric 
matrices introduced in our previous work [20]. For more details on RMT applications 
to the QCD-like Dirac operator spectrum we refer to [21J. 

Many more non-Hermitian RMTs than just the three Ginibre ensembles and their 
chiral (or Wishart /Laguerre) counterparts exist [22] and these are mostly unsolved 
to date. Very recently another two-matrix model generalisation of the real Ginibre 
ensembles was introduced and solved in |23j. There the eigenvalue correlations of the 
ratio of two quadratic matrices are sought, whereas here we deal with the product of two 
rectangular matrices. Whilst the former case leads to a Cauchy-type weight function, in 
our model we will obtain a weight of Bessel-.fr functions for the eigenvalues. We hope 
that given the plethora of RMT applications, our particular model will find applications 
beyond the field theory that it has been designed for. 

The approach of solving our model is based on the variational method detailed in 
[HUD]. It follows its two main ideas: first to compute the joint probability distribution 
function (jpdf) for general N by reducing it to 2 x 2 and lxl blocks. Because we are 
considering rectangular matrices this is a priori not guaranteed to work. Second, we 
use the variational method [7J [10] in combining all density correlations with n points 
(being real, purely imaginary or complex conjugates) into a single Pfaffian form. This 
reduces the computation to determining its main building block, an anti-symmetric 
kernel. Whilst it can be deduced from the spectral 1-point density - which was known 
for the real Ginibre ensemble [21] - we here exploit an idea from our previous publication 
|20j. There the kernel was determined by computing the expectation value of two 
characteristic polynomials using Grassmannians. The same relation between kernel and 
characteristic polynomials is known to hold for the symmetry classes with complex [25] 
or quaternion real matrices [19] . in fact for any class of non-Gaussian weight functions. 

As a new result we can express our kernel for real asymmetric matrices in terms 
of the kernel for complex non-Hermitian matrices. Such a relation might have been 
expected to exist as it is known for Hermitian RMT [261 [27] . 

Other methods that have been applied successfully to the real Ginibre ensemble such 
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as the supersymmetric method skew-orthogonal polynomials [5] or probabilistic 
methods [9] are very likely to be extendible to our two-matrix model as well. 

The paper is organised as follows. In Section [2] we summarise our main statements: 
the definition of the matrix model, its jpdf in terms of the real, imaginary and complex 
conjugate eigenvalue pairs, and the solution for all density correlation functions as a 
Pfaffian of a matrix-valued kernel. Examples are given for the simplest spectral densities 
at finite N and in the microscopic large- N limits for strong and for weak non-Hermiticity 
at the origin. These findings are then detailed in Section [3] on the jpdf, where we 
separately treat N = 1, 2 and general N. The spectral density correlations and their 
finite- and large- iV results are derived and illustrated in Section HI Our conclusions are 
presented in Section [5j Some technical details on the computation of the Jacobian are 



collected in Appendix A 



2. Summary of results 

2.1. The model 

The chiral Gaussian ensemble of real asymmetric matrices as introduced by the authors 
[20J is given by a two- matrix model of rectangular matrices P and Q of sizes A" x (N + u) 
with real elements, without further symmetry restriction. The partition function 
normalised to unity is given by 

= J / dP dQ exp - 

where we integrate over all the independent, normally distributed matrix elements of P 
and Q. We are interested in the eigenvalues of the matrix T> of size 2 A" + v squared 



2 -Tr(PP T + QQ 1 



(2.1) 





I 




H 




1 



»=\ ~t ° . P Y Q i-i ; T " i • (2.2) 



P L - nQ 

Here fi G (0, 1] is the non-Hermiticity parameter, interpolating between the chGOE 
(lim/i — > 0) and maximal non-Hermiticity (fi = 1). The analogous chiral Gaussian 
two-matrix models with unitary and symplectic symmetry were introduced in [THl fTU] 
respectively. 

In applications to field theory, T> corresponds to the chiral Dirac operator, and fi to 
the chemical potential Typically, Nf extra determinants of the type det[D + ml2N+u] 
are inserted into the partition function eq. (12.11) . where m is the quark mass, but we will 
restrict ourselves in this paper to the case Nf = 0; this is referred to as the quenched 
case. 

For later convenience we give an equivalent form of eq. (12.1 p , by changing variables 

from 

P=\(A + B), Q = ±-(A-B), (2.3) 

| The Euclidian Dirac operator in field theory is actually anti-Hermitian for n = 0, but we will not 
use this convention here. 
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to matrices A and B defined in eq. (12.21) : 

/ -. \ N{N+u) r r 
Z=[—) / dA dB e -h + ^(AA- + BB-) + ^TT(AB- + BA-) ^ {2A) 

with r]± = ^ 2 . (2.5) 
The two /z-dependent combinations rj± will be used throughout the paper. 



2.2. Eigenvalue representation 

The eigenvalues A of the Dirac matrix D are determined from the following equation^: 

N 

= det[AI 2N+u - V] = A" det[A 2 I N - AB T ] = A u JJ(A 2 - A 2 ) . (2.6) 

i=i 

For this reason we will first compute the eigenvalue distribution of the N x N Wishart- 
type combination of matrices C = AB T . C has real elements, and therefore its 
eigenvalues A 2 are real, or else come in complex conjugate pairs. The matrix T> itself 
has the following solutions: v zero-eigenvalues A = 0, and 2A^ eigenvalues coming in 
pairs A = ±Aj. Consequently the non-zero eigenvalues of T> fall into three categories: 

(i) for A 2 > 0: real pairs A = ±Aj e K 

(ii) for A 2 < 0: purely imaginary pairs A = ±A 3 - G iM. 

(iii) for pairs A 2 , A* 2 : quadruplets A = ±A i? ±A* e C \ {R U iR} . 

This leads to an accumulation of eigenvalues on both the real and the imaginary axes 
as already pointed out in [2D] . The same phenomenon has been observed numerically in 
a one-matrix model [Hj based on the proposal [13] (obtained from eq. (12.11) by choosing 
Q ~ I). This is in contrast to the real Ginibre model where eigenvalues accumulate 
only on the real axis (sec e.g. [29|). 

The joint probability distribution (jpdf ) for the matrix C is obtained from eq. (12.41) 
by inserting a matrix delta function; using the cyclic property of the trace we then have 

P(C) ~ exp[r?_TrC] / dA f dBexp \——Ti(AA T + BB T )] 5{C - AB T ) . (2.7) 

lmN(N+v) 1m,N(N+v) L 2 



Zk = %k + iyk = A 2 with d 2 z k = dx k dy k is 



As shown in Section [31 our final result for the jpdf of V in terms of squared variables 

Z I d 2 Zl ... I d 2 z N P N (z u ...,z N ) (2.8) 
Jc 

N P N [N/2] 



C 



J*n _iV L / J > h 

jv[] / d 2 z k w{z k ) \\{zi - Zj) ^ (T[(-2i)5(x2i-i - X2i)6(y 2 i~i + y2i)®{y2i-i) 

k=l^ C i<j n=0 1=1 

X Q(X! > X 3 > ... > X 2 n-l)Q(x 2 n+l > %2n+2 > ■■■ > X N )5(y 2 n+l)-S{yN)j- (2.9) 



This follows from eq. (|2.2p using the standard relation dot ( ° ^ ) = det(d) det(a — bd 1 c) for 



c 



square matrices a and d, with d invertiblc. 
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The integration measure d 2 Zk extends over the complex plane for each of the z^. The 
normalisation constant cjy will be given in eq. (13.461) . In eq. (12.91) we sum over all 
distinct possibilities for N eigenvalues to come in n > complex conjugate pairs, with 
the remaining N — 2n > eigenvalues being real. For n = in eq. (12.91) the product in 
the first line is simply unity. 

Specifically, the jpdf is only non-zero when the eigenvalues appear in the following 
order: the n complex eigenvalue pairs must be placed first, ordered with respect to 
decreasing real partg[jj] 

(3mA 2 > 0), (KeA 2 = KeA 2 , 3mA 2 = -3mA?), (KeA 2 < 3fteA 2 , 3mA 2 > 0), . . . , 

(3?eAL_! < 3?eAi n _ 2 ,3mAL-i > 0), (KeA 2 „ = HeA^-i, 3mA 2 n = -SmA^^), (2.10) 

and the N — 2n > real eigenvalues follow, and are also ordered with respect to 
decreasing real parts: 

AL+i > AL+2 > ••• > A^ . (2.11) 
The function g(z) inside the weight function 

w{z) = \z\ u/2 exp[rj-z] g(z) (2.12) 
depends on whether z (= A 2 ) is real or complex: 

g{z) = 2Ke( V+ \z\) , for z GM , (2.13) 

roc t. 

[g{z)f = [g{z*)f = 2 / e^^-v 2 )-** K* (2 V 2 + t(x 2 + y 2 )) edc(2ri + y/i\y\), 

Jo 1 

for z = x + iy G C. 



(2.14) 



Because the complex eigenvalues come in pairs we will always get two factors g(z) for 
each pair (note the square on the left hand side of the definition eq. (I2.14p ). The limit 
y — > of a single g(z) in eq. (12.14ft is smooth, leading to eq. (12. 13ft . 

The essential idea in the derivation of the jpdf detailed in Section [3] is to reduce 
the calculation of the jpdf for D with general N down to 2 x 2 and lxl blocks, which 
can be handled in terms of the N = 2 and N = 1 problems which we solve explicitly. 



2. 3. Density correlation functions for finite N 

We follow the method of using a generating functional for all eigenvalue density 
correlation functions introduced in [3, [10] . Because we essentially follow 10] we can 
be brief here. We enlarge the definition of the partition function eq. (I2.9P by introducing 
sources /(A|)0 

Z[f\ = / d 2 z x ... I d 2 z N P N (z 1} z N )f( Zl ) . . . f(z N ) . (2.15) 

j Unlike for real eigenvalues two complex eigenvalues that are not complex conjugates can have the 
same real part, without the jpdf vanishing. Although being of measure zero we can fix this ambiguity 
by ordering with respect to decreasing absolute imaginary part. 

These will symmetrise the ordered eigenvalues when differentiating. 
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For pairwise distinct arguments Z\ 7^ z 2 7^ . . . 7^ z n the n-point density correlation 
functions are then generated in terms of functional derivatives with respect to the 

5 2 {z-z'): 

5 5 



sources, leading to insertions of delta functions S j, 



-^n(^l) ■ ■ ■ ; ^n) 



Sf(z) 
Sf(z') 

-z[f\ 



(2.16) 



5f( Zl ) ■■ m sf{z n y 

In doing so, each n-point function contains a sum of different contributions, splitting n 
into all possible combinations of real eigenvalues and complex eigenvalue pairs. Both 
the generating functional Z[f] and the n-point density R n can be written as Pfaffians 

mm- 



Z[f] = c N Pf 
]{, n (21 . . . . , z n ) 



d 2 Zl / d 2 z 2 f(z 1 )f(z 2 )T(z 1 ,z 2 )zl- L 4 



i-l J -I 



(2.17) 



l<i,j<N 



(N-n)\J c 
Pf 



n+l 



d 2 z N P N ( Zl , 



zn) 



K, N (zi,Zj) -G N ( Zi , Zj ) 
G N (zj, Zi ) -W N (zi,Zj) 



(2.18) 



l<i ,j<n 

In the latter case, one has to compute the Pfaffian of the ordinary, 2n x 2n matrix 
composed of the matrix of quaternions inside the square bracket. We have restricted 
ourselves to even N for simplicity. The case of odd N can be treated along the lines of 
[TP] (or [TTJ for an alternative formulation). We have introduced the following functions 
of two complex variables Zj = Xj +iyj, j = 1, 2: 

w(z l )w(z 2 )(2i5 2 (zi - zl)sgn{y x ) + %i)% 2 )sgn(x 2 - ari)V (2.19) 

V- 



T(z x ,z 2 ) 
K,n{z\,z 2 ) 



8tt(4uV)^ +1 



0=0 



V- 



8vr(4uV)^ +1 







Z 2 T\ Z\— — 

OZ 2 OZ\ 



ju ( *2 

J+1 V4uV 
d z 2 - 



Z 2, 



(2.20) 



Z\ 



L 



z 2 



A^ 2 r]_ 



G N (zi, z 2 ) 
WnUi, z 2 ) 



(j + v)\ 3 V4u 2 ?7„ 
j— « 

/ d 2 zK, N (z 1 ,z)F(z,z 2 ) , 
Jc 

J r (zi,z 2 )+ / d 2 z / d 2 z'J 7 (zi, z)JCn(z, z')T{z\ z%) . 
Jc Jc 



(2.21) 

(2.22) 
(2.23) 



The kernel JCn{zi, z 2 ) is the building block of all correlations, and is actually defined in 
terms of T{z\,z 2 ) (see Section H] for further details). However, we will give a different 
proof of the precise form of eq. (I2.20p which does not rely on direct evaluation of the 
definition; rather, it will be calculated from the expectation value of two characteristic 
polynomials (see also [2U]). Note that in eq. (12.211) we have expressed the kernel of our 
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real two- matrix model {(3 = 1) as a derivative of the kernel of the complex two-matrix 
model [16] (/3 = 2). We found that a similar relation holds relating the kernel of the 
(3 = 1 real Ginibre ensemble [HI [20] to the one for the f3 = 2 complex Ginibre ensemble 
[30] , containing Hermite polynomials in the complex plane. 

There is an integration theorem in analogy to other matrix models with complex 
eigenvalues [3 EH [26] : 



d 2 z n Pi 



(iV-n + l)Pf 



^■N{Zi,Zj) —GN{Zi,Zj) 

G N (zj,Zi) -W N (zi,Zj) 

J l<i,j<n 

fC N (zi,Zj) ~G N (zi,Zj) 
G N {zj,Zi) -W N (zi,z s ) 



(2.24) 



l<i ,j<n— 1 

This is the major advantage of working with an n-point correlation function as defined 
in eq. fl 2 . 1 8 [) which contains all possible contributions of real eigenvalues and complex 
conjugate eigenvalue pairs. Had we studied instead a particular n-point function with a 
fixed and given number of real and complex conjugate eigenvalues, then we would have 
found that such a simple integration theorem does not exist [5]. 

Let us spell out one example explicitly, the spectral density, which we will give for 
even iV only. Details can be found in Section H] including figures, where we follow the 
method of [3, \10\ . From eq. (I2.18j) we have 



Ri{zi) 



d 2 z jc n ( Zi , z)f(z, Zl ) = r^{ Zi ) + %iX0n; 



(2.25) 



Inserting the appropriate weight from eqs. (12. 13[) and (12.141) we obtain with z = x + iy 

R^{z)= -2i\z\ u e^- x sgn{y)K, N {z 1 z*) (2.26) 



x 2 



dt 



exp 



2 V 2 + t(x 2 - y 2 ) - 1] K» (2r, 2 + t{x 2 + y 2 ))erfc ^v^l 



Rf(x) 



dx'sgn(x 



f 2e v-(^')2K»( v+ \x\)2K»( v+ \x'\)lC 



N {x,x 



(2.27) 



Eqs. (I2.25P - (I2.27P are valid for even iV only. In the final step we can change 
from squared variables z = x + iy = A 2 to Dirac eigenvalues A, using the simple 
transformations 



RlDirac( Z ) 

Rf 



Dirac 



A\z\ 2 Rf(z 2 ) 
2\x\Rf(x 2 ). 



(2.28) 



Note that the latter describes the density both of real Dirac eigenvalues, for Rf£> irac (x) 
with x G M, and of purely imaginary Dirac eigenvalues, for Rfuirac( x ) with i G i8. 
Because Rf(x 2 ) ^ Rf(—x 2 ) this is not the same function, see e.g. Fig. [2] in Section [H 



2.4- The large-N limit 

In order to take the large- N limit in principle one first has to rescale all eigenvalues in 
eq. (12.91) A/. — * \/N A&, which is equivalent to giving all the matrix elements in eq. (12.11) 
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Gaussian weights exp[— (N/2)TrP T P] for P, and similarly for Q. In this parametrisation 
the macroscopic spectral density will have a compact support in the large- N limit, given 
by a circle for \i = 1, and an ellipse for < // < 1. We will not discuss this macroscopic 
limit in detail but will focus on the local correlations, i.e. the microscopic large- N limit. 
Here one has to distinguish between strong and weak non-Hermiticity [32] , each of which 
involves a second rescaling. 

We first give the strong non-Hermiticity limit after the first rescaling, defined by 
keeping \i fixed and only rescaling the eigenvalues according to 

A 2 . (2.29) 



lim NA 2 



This scaling actually cancels the first scaling. Because of this the scaling limit is also 
true away from the origin. 

We only give the microscopic kernel here, to be inserted into eqs. (I2.26P and (I2.27p . 
For fj, = 1 the kernel simplifies to monic powers, see eq. (33) in [20] . Its large- N limit is 
easily seen to lead to a modified Bessel function 



Zl,Z2 



,Zi - z 2 ) 
Q4n2» 



E 

j=0 



/£l£2\: 

j\(j + is)\\ 4 J 



Z\ ~ ^2] 

64vr 



[Z1Z2 



Z1Z2 



(2.30) 



Because /1 is not scaled here, in contrast to the weak limit below, the insertion into eqs. 
(I2.26P and (12.271) is relatively straightforward, apart from the phase for negative real 
eigenvalues. The case for general /1 which can also be obtained by simply rescaling the 
arguments in eq. (12.301) is discussed in Section 0] where we also show plots. 

The weak non-Hermiticity limit at the origin (again after the initial rescaling above) 
is defined by scaling both the squared Dirac eigenvalues A 2 and the chemical potential 
/i 2 with 2A^, corresponding to the volume in field theory: 

2N{i 2 = a 2 , lim (2N) 2 A 2 = A 2 . (2.31) 



lim 

N—>oo,fi- 



In this limit the macroscopic density is projected back onto the real axis, with probability 
density given by the semi-circle as for fi = (when our model is Hermitian), whilst 
microscopically the eigenvalues still extend into the complex plane. 

The limiting microscopic kernel, as a function of squared variables z = A 2 , can be 
expressed in terms of our completely unsealed finite- kernel as 



1C W \ 



z l, Z2, 



= lim 



( ZlZ2 



(4A0 2 V(4A0 2 
•l 



v/2 



Y~ 1 Z\ Z 2 



a 



*2N 



ds s 2 * 



-2a 2 s 2 



lrr . , {y/zl Ju+i(sy/zT)Ju(sy/z2~) - (zt <-> z 2 )} (2.32) 

which is to be inserted into the definition of the density eq. ( 12. 25ft . From eq. ( 12.261) we 
obtain for the microscopic density of complex eigenvalues 



— 2i sgn(y) exp 



x 

I4c7 2 



K 



w , 



z, z 



x 2 



dt 

T exp 



t{x 2 - y 2 ) 
32a 4 



1 

4tJ 



2 



t(x 2 + y 2 ) 
32a 4 



erfc (i! 2 



y\ 



(2.33) 
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The real eigenvalue density is more subtle; we again refer to Section H] for a more detailed 
discussion of the weak limit, including figures. 

3. Calculation of the joint probability distribution function 

In this section we compute the joint probability density function (jpdf) as stated in 
eq. f)2.9p for the squared non-zero eigenvalues of T>. For pedagogical reasons we first 
compute the jpdf separately for N — 1 and 2 in Sections 13.11 and 13.21 respectively. 
This is because we will need these results when treating the general N case in Section 
13.31 as these sub-blocks will appear in the computation of the general Jacobian. Some 



technical details will be deferred to Appendix A The cases with N — 1, 2 will make the 



parametrisation and residual symmetries more transparent for later. 
3.1. The N = 1 case 

In this simplest case our matrices P and Q, or after changing variables A and B in eq. 
( 12. 3p . are of size 1 x (1 + v) and are thus given by vectors a and b, each of length u + 1. 
The eigenvalue equation (12.61) for T> becomes 

= A"(A 2 -a-b) , (3.1) 

and thus we only have a single non-zero (and real) eigenvalue A 2 to determine. Its (j)pdf 
is given by 

P(A 2 ) = - 1 —— [ da / alh exp[??_a • b] exp [-5+ (| a | 2 + |b| 2 )l S(A 2 - a • b) . (3.2) 

(47r^r +1 Jw+i Jwv+i ' L 2 J 

We simplify this expression in two steps. First, without loss of generality we may choose 
the direction of b as the first basis vector for a in Cartesian coordinates. This leads 
to a decoupling of the remaining components a 2 , . . . ,a u +i, and the integral now only 
depends on b through its modulus b = |b|. Second, we choose polar coordinates for the 
vector b, leading to the Jacobian b u , and, on symmetrically extending the integral over 
b to — oo, we obtain 

1 9--n-( v +i)/ 2 /q-7t\ ^/ 2 i r°° r°° 

^ - (4^4 m (£) ^i/^/^e-KC*^-^ 

e"- A2 / dbe-^+^lbl"- 1 . (3.3) 



23^/2+2^+1^/2^ r 

The first new pre-factor comes from the surface area of the unit ^-sphere 
_ 27r(-+D/ 2 VO(u+l) 

(with VO(v) being the volume of the orthogonal group) through the angular integration 
over b, the final pre-factor from the Gaussian integrations over the decoupled 
components of a. It is important to note that the first line of eq. (13.31) looks like a 
reduction to the v = case, apart from the extra factor \b\ u from the Jacobian. We will 
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use the same strategy for iV = 2 in the following subsection. In the next step we change 
variables e* = 6 2 /|A 2 | to arrive at 



i r°° 

P(A 2 ) = |AIV- A2 2 / coshfz/^e-^ 1 ^ 1 ^ 11 *^ 

2 3-/2+2^+l^/ 2 v ^ r (£±l) 1 J 

-|A|V>- a2 2A%(?7+|A 2 |) . (3.5) 



Here we have used a particular representation eq. 9.6.24 in [33] of the i^-Bessel function. 
It directly gives cn=i times the weight function in eqs. (I2.12j) and (12. 13f) when changing 
variables A 2 — > z, where 

c n=i = ; ~p? TT- (3-6) 

This is consistent with eqs. ( 12. 9p . (12 . 12[) and (I2.13p . and ends our calculation for N — 1. 

As a remark, eq. (13. 5p can be derived in various different ways, including by Fourier 
transformation. It is known that if a and b are independent random variables with 
normal distributions, then the product c = ab has distribution function P ~ K . 
Consequently, the sum of v + 1 (independent) such variables q has a distribution given 
by the convolution of v + 1 functions Kq. Fourier transformation F turns this into an 
ordinary product, and so we obtain P v = F~ l {(F(P )) U+1 } ~ K», i.e. on performing 
the integrals we again reach eq. (13.51) . 



3.2. The N = 2 case 

Our matrices A and B are now given by two row vectors a J=1) 2 and b J=1) 2 each of length 
z/ + 2: 

A=( ai ) , B=(* 1 ) , C = AB T =( ai ^ 1 ai ^ 2 V (3.7) 
y a 2 J y b 2 J y a 2 ■ bi a 2 • b 2 J 

The eigenvalue equation (12.61) 

= A u det[A 2 1 2 - AB T ] (3.8) 

has two solutions A 2 2 which may (i) both be real or (ii) form a complex conjugate pair. 
We will have to distinguish these two cases below. 

For the first step we reduce the calculation of P{C) in eq. ( 12.7ft to a matrix integral 
of 2 x 2 matrices A' and B' times a Jacobian, as we did for = 1 in the first line of 
eq. ( 13. 3p . The resulting Jacobian here will be ~ | det B'\ v . Our aim is to rotate both b x 
and b 2 into the x?/-plane of the coordinates for the a.j, as then 

v+l 

a * ' b i = aikb i k = + ai2& J' 2 > ( 3 - 9 ) 

k=l 

and the remaining components of the a, decouple. Here we use primed vectors and 
coordinates to denote the quantities after rotation. 
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For v = 1 the corresponding Jacobian is obtained as follows. Rotating the 3D vector 
into 2D by h[ = O hi gives rise to a factor ~ |bi|. Alternatively it can be computed 
by comparing the initial and final 'volumes' (generalised surface areas, in fact), yielding 
^l? 1 ^ = 2|bx|, where S n is given by eq. (13.41) . The remaining rotation around h 1 rotates 



b 2 into the xy-plane as well, giving 



5i|b 



2 1 smt 



So 



7r|b 2 |sin#, in which 9 G [0, tt] is the 



angle between b x and b 2 . The final Jacobian reads 
| hi | | b 2 1 sin 9 = \h'i \ \h' 2 \ sin 9 - 



b'n b' 12 
b' h' 



det B'\ 



(3.10) 



where the last equality easily follows by parametrising h[ and b' 2 in 2D, and the factor 
|£ corresponds to the angular integration over b. 

For v > 1 we can thus successively repeat these steps by projecting onto one 
dimension lower until we reach the xy-pl&ne. The volume factors will telescope out and 
we arrive at 



Si/Siz+i 
SqSi 



( | h[ 1 1 b 2 1 sin l 



SqSi 



det B' 



(3.11) 



We thus have reduced eq. (12.71) for N = 2 from 2 x (2 + v) down to 2 x 2 matrices: 



P(C) 



c N=2 y+ „_Trc 



2tt 3 



(3.12) 



x 



dA' / dB' exp 



/ AlT 



-- V+ Tr(A'A 



B'B 



5(C-A'B' T ) \detB'\ 



where ctv =2 is defined in anticipation of the final result as 

° N=2 ~ 8tt (2vr)-(2 / i)^;+i Sq Si ' [6A6) 

The Gaussian integrals over the decoupled components of the two vectors ajk for k > 2 
have been evaluated, using that TtAA t = \sli\ 2 + |a 2 | 2 . The 2x2 matrix A' can now 



be integrated out, by formally changing variables A' 
I det B 



A'B with Jacobian 



/|-2. 



P(C) = ^±e^ c [ dB' exp 



■~ V+ Tr (CC T (B'B' T )^ 



B'B' T ) 



det B 



/ 1 iv— 2 



(3.14) 

Note the similarity with eq. (13. 3p . 

In a second step we perform the integral J dB'. Because CC T is a symmetric, 
positive definite matrix we can diagonalise it with an orthogonal transformation O 



T (CC T )0 




, A li2 >0 



(3.15) 



Using the invariance of dB' we can change variables B' — > OB' with (B'B' T ) 1 — > 
0(B' B' T )~~ 1 T . We thus replace CC T by its diagonalised form eq. (13.151) in the exponent 
in eq. (13331) : 

Tr (B'B' T + CC T (B'B IT )- 1 ) = a 2 + b 2 + c 2 + d 2 + h- 2 ((c 2 + d 2 )\i + (a 2 + 6 2 )A 2 ), (3.16) 
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(3.17) 



We now introduce h as an independent variable in eq. (I3.14p by inserting a delta-function 
constraint in its integral representation: 



P(C) 



c N=2 y+ „Trc 
4vr 4 

x exp 



dh\h\ 



u-2 



dB' 



du! e -Mh—(ad-bc)) 



(3.18) 



— oo 
-2//J2 , j2> 



-^[a 2 + b 2 + c 2 + d 2 + h~ 2 {{c 2 + d 2 )A a + (a 2 + 6 2 )A 2 )] 



Cjy =2 77+ „Trc 
4vr 4 



^-2 



(27T) 



<L; e 

w 2 + 77 2 (l + Ai/ ^2 )(1 + A2//l 2 ) 



where we performed the Gaussian integrals successively in pairs a, d and b, c, and 
switched to positive h. The denominator in the third line can be rewritten at the 
cost of an additional integral, | = J °° dte~ at for a > 0, and after changing variables 
uj — > r = cu/i we have 



P(C) 



2c ^v=2 ?7+ « Tr 



c 



7T 



2cat = 2 ^+ 



7T 



3/2 



(AiA 2 )"/ 4 e* 



dh h 



Tic 



u-l 



dre 



dt 

7! exp 



( fc e -(T 2 W+(h 2 +(M+\2)+\l\2/h 2 ))t 



-^+(Ai + A 2 )t - - 



4f 



K f (2r ? 2 v / A7A 2 "t). 

(3.19) 



We performed the Gaussian integration over r first, and then employed the following 
integral from [M] eq. 8.432.6 

-t-z 2 /4t 

(3.20) 



1 , £ 

2 V2 



dt 



to do the h- integration after another change of variables, arriving at our result eq. (I3.19p . 

As the final step we need to express P(C) in terms of the eigenvalues of C (i.e. A 2 
and A 2 ,), rather than those of CC T . This will also involve a Jacobian to be computed 
later. In eq. (I3.19P we need 

TrC = A 2 + A 2 , AiA 2 = det[CC T ] = k\k\ , (3.21) 

which are trivial. Only the combination Ai + A2 = Tr(CC T ) requires more calculation. 
In general, we can orthogonally transform any real matrix C to the form 



C 



sin 9 — cos 9 
cos 9 sin 9 




sin 9 cos 9 
- cos 9 sin 9 



(3.22) 



where e\, €2, s and the rotation parameter 9 are all real. The matrix parameters {e\, e 2 , s} 
and eigenvalues {A 2 , A 2 ,} follow from 







ei-A 2 



e 2 -A 2 



A 



{e x + e 2 )A 2 + (eie 2 + s 2 



(3.23) 
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with solution 



A 2 



ei + e 2 



± 



4 



This can be inverted for f° r later use to 

± 



e l,2 



A 2 + A 2 



(A? 



A 2 .) 2 



+ s 2 



(3.24) 



(3.25) 



2 V 4 

This immediately leads to 

A a + A 2 = Ti(CC T ) = e\ + e\ + 2s 2 = A* + A A 2 + 4s 2 . (3.26) 

On inserting eqs. (I3.2ip and (13.261) into eq. (13 . 1 9[) we can express P(C) in terms of 
A 2 2 . However, we have changed variables twice to arrive here, first from the matrix 
elements Cjj of C to {ei, €2, s, 9}, and second from {ei, e 2 } to {A 2 , A 2 ,}. The corresponding 
Jacobians we must multiply by are given by 
d{c n ,c 12 , c 2 i,c 22 } 



Ji 
Ji 



d{e 1 ,e 2 ,s,9} 
d{e 1 ,e 2 } 



2 |ei - e 2 | cos 2 (2#) 



9{A 2 ,A 2 } 





A? 


-A3 




ei 


- e 2 



(3.27) 



3.2.1. Distinction of real and complex eigenvalues 

In order to give the jpdf from eq. (13.191) for variables A 2 2 alone we have to integrate 
over the remaining real variables s and 9. Here we have to distinguish between the case 
of two real eigenvalues, and that of a complex conjugate pair. Starting from a real 
matrix C in eq. (13.221) all new variables, in particular ei j2 , are real. In view of eq. (13.251) 
there are two possibilities if the radicand is to remain positive: 

(i) A 2 , A 2 , are both real =>- \(A\ — A 2 .) 2 + s 2 > which is always satisfied, or 



(ii) A 2 , A 2 , are complex conjugates: A 2 — A 2 6 i 
We thus obtain for the jpdf 

dAldA 2 2 P(AlA 2 2 ) = 2CN=2V+ 



s 2 >-m 



A 2 ) 2 > 0. 



x 



3/2 



ds2 



dA\dA 2 2 e^ +K V\AlAl\ v l 
00 dt 



2tt 



Vt 



cxp 



nUA\ + At + ^)t- 



cft?|A 2 -A 2 |cos 2 (2fl) 
1 



At 



K^2rf + \AlAl\t), 
(3.28) 

where s 2 min = max{0, —(A 2 — A 2 ) 2 /4}. The ^-integral is trivial. 

For two real eigenvalues having s m i n = the integral over s can be performed, 
leading to the following simplification for the remaining t-integral: 



(3.29) 



2 f j e-^ (A ' +A ^"^ f (2 V 2 + |A 2 A 2 |t) = 2K~ (^ + |A 2 |) 2^(^ + |A 2 |) . 

Here we have used [34J eq. 6.653.2 after changing variables t — » u = ^. When ordering 
the two eigenvalues as A 2 > A 2 , the jpdf eq. ( 13.281) can thus be written as 

P(A 2 , A 2 ) = c N=2 (A 2 - A 2 ) J] |A 2 r/V" A ? 2K»( V+ \A 2 \), A 2 , 2 e R , (3.30) 

i=i,2 
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as was claimed in eq. (12.91) in conjunction with eq. (12.131) for N = 2. 

For two complex conjugate eigenvalues the s-integral leads to the complementary 
error function, without further simplification. Here we order < SmAj = — SmA 2 , as 
in eq. (12.101) to obtain the following real, positive distribution 

P(A 2 , A 2 2 ) dA\dA\ = c N=2 dA\dA 2 2 (A 2 - A 2 2 )\A 2 ^/ 2 \A 2 2 \^ 2 e^ A ' +A ^ (3.31) 

POD 7. 

x2 / - e -^^ +A ^-^ir f (2r / 2A2A^)erfc(r ?+v / t|A2-A2|), A\ 2 e C . 
Jo t 

Because the integral depends only on the modulus |SmA 2 | = l^mA^I we can define its 
square root to be the weight of each eigenvalue A 2 , as in eq. (12.141) . The limit l^mA 2 ] — > 
smoothly reduces the integral in eq. (I3.3ip to eq. (I3.29p . using that erfc(0) = 1. 

Combining the real and complex cases, and still assuming the eigenvalue ordering, 
we can write the jpdf most generally as follows 

d 2 z 1 d 2 z 2 P(zi, z 2 ) = c N=2 d 2 z l d 2 z 2 (zi - z 2 )w(z 1 )w(z 2 ) (3.32) 
x (s(y 1 )5(y 2 )Q(x 1 - x 2 ) - 2id 2 { Zl - z* 2 )Q{ yi )) 

where we have switched variables A 2 = z = x + iy, including for the differentials, 
dA 2 dA* 2 = (dx + idy)(dx — idy) = —2idxdy = —2id 2 z. This is just the jpdf in eq. ( 12.91) 
including the weight eq. (12.121) . and this completes our computation for the N = 2 case. 



3.3. General structure for arbitrary N 

In this subsection we will compute the jpdf for any - both even and odd - given in 
terms of the eigenvalues A 2 of the Wishart matrix C = AB T . We start from eq. (12.71) 
which we repeat here for convenience, 

P(C') ~ ex]>[//_IrCj / <1A I dB exp -^Tr(AA T + BB T ) 8{C - AB T ) . (3.33) 



N+v) JrjlN(N+v) 



We will eventually use the results from the previous two subsections as building blocks. 

Instead of using a generalised Schur (or QZ) decomposition involving unitary 
matrices to bring A and B to upper triangular form, we will restrict ourselves here 
to orthogonal transformations. The best we can achieve in this way is a so-called 
almost (or quasi) upper triangular (AUT) form for one of the matrices, and an upper 
triangular form for the second (which is also AUT) . An AUT matrix is composed of a 
block diagonal matrix, having non- vanishing 2x2 blocks along the diagonal for even N, 
and an additional lxl block at the end if is odd. The remaining non-zero elements 
of an AUT matrix all lie above the block diagonal. 

The precise transformation that we make is: 

A = A (A A + A A )0%, B T = B (A B + A B )0 T A . (3.34) 

Here A^ and A# are block diagonal, and and A# are zero except in elements strictly 
above the block diagonal. A a + A^ and A B + A B are hence AUT. Note that Oa is of 
size N x N, and B is of size (N + v) x (N + v). A a and A^ are each the same size as 
A itself (i.e. rectangular), and similarly A# and A B are the same size as B T . 
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Matrix 


Degrees of freedom 


Matrix 


Degrees 
Even N 


of freedom 
Odd N 


A 


iV(iV + v) 


Aa 


2N 


2N-1 


B 


N(N + v) 


Ab 


2N 


2N- 1 






A A 




Ar2 + 1 N + vN 






A B 


2 


N \ +l N 

2 






o A 


2 


JV2 + 1 N 






B 




N \ +1 N + vN 



Table 1. Counting dof before and after change of variables. 



To make the transformation eq. fj3.34j) unique, having the same number of degrees 
of freedom (dof) on the left- and right-hand sides, we restrict O A and Ob as follows: 

A G 0{N)/0{2) N/2 , B G 0{N + u)/0{2) N/2 0{u) for even N , 

A G 0(N)/0(2)^ N - 1 '>/ 2 , Ob G 0(N + u)/0{2)^- 1 ^ 2 0{y) for odd N. (3.35) 

The residual symmetries leading to these cosets are rotations within each block on the 
diagonal, and, loosely speaking, the extra factor 0{y) can be thought of as originating 
from the reduction of B from rectangular to square form, as in the N — 1 and 2 cases 
earlier. The precise counting of dof is given in Table [H matching the sum of dof of A 
and B for all N. 

Under this orthogonal transformation, the integrand in eq. ( 12. 7ft changes as follows: 



exp 



-yTr (AA T + BB T ) 



exp 



-^Tr(A A A^ + A B A T B + A A A T A + A B A T B ) . (3.36) 



remains unchanged, with the relation 
es yet to be determined. 



The pre-factor exp [^_TrC] = exp 7/_ YliLi Af 
between the eigenvalues A 2 and the new variab 

The transformation eq. (13.341) leads to the following differentials 

(<L4)y = (o A [OldO A (A A + A A ) - (A A + A A )0 T B dO B + dA A + dA A ]0 T B 
(dB T ) tJ = (0 B {0 T B d0 B {AB + Ab) - (A B + A B )0 T A dO A + rfA B + dA B }0 T A 

Here we have used the fact that for orthogonal transformations (with T = I) the 
differential T dO is anti- symmetric. This remains true for our special choice of cosets. 
When considering the invariant line element Tr(cL4 dA T + dB dB T ) the rotations outside 
the square brackets in eq. ( 13. 37ft can be dropped. 



(3.37) 
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We will now compute the Jacobian for the change of variables from {dA, dB T } to 
{dAji, dAs, dAA, dAs, A dOA, O^dOs}- Here we use the differentials for convenience as 
for the orthogonal matrices only these constitute independent variables (see also [31] for 
a similar discussion for the real Ginibre ensemble). In particular, counting dof T A dOA is 
an anti-symmetric matrix with zeros not just on the diagonal but on the block diagonal. 
Similarly Oj^dOs is an anti-symmetric matrix with zeros on the block diagonal for the 
first N elements, and a zero-block of size v x v on the remaining part of the diagonal. 

In the differential eq. (13.371) the variables {dAA, dAs, dAA, dAs} are already 
diagonal, with 



d(dA) tJ djdA^ 
d(dA A ) pq '" ' d(dA " " ' 



$ip$jq i <v J A \ — $ip$jq > (3.38) 
A)pq 



and similarly for dB T . This contributes a unity matrix block to the Jacobi matrix 
J , when considering the corresponding elements of {dA, dB T } on and above the block 
diagonal. 

The non-trivial contribution from the Jacobian therefore originates from 
differentiating the remaining elements of {dA, dB T } below the block diagonal with 
respect to the independent variables of {0 T A dO 'a, O^dO b} , where we also choose the 
lower block diagonal elements for convenience. When appropriately ordering the 



elements of the Jacobi matrix (see Appendix A and also a similar discussion in [31] ) 



the contributions proportional to and A# in eq. (I3.37P will drop out, being part of 
a lower triangular sub-matrix in J . 

For the sake of argument we restrict ourselves in this section to the case of A^ and 
Ab being diagonal. The more general (and typical) case in which A^ and Ab contain 



2x2 blocks is treated in Appendix A 



Arranging the remaining matrix elements below the block diagonal of dA and of 
the square part of dB T into pairs this leads to a 2 x 2 block diagonal Jacobi sub-matrix 
with elements 

d{dA)jj d{dB) i:j 



Hpf d{0 1 A dO A ) l] diO^dOAhj 
<l<1 I d{dA)jj d{dB)jj 



det 



\A)-Al\ , (3.39) 



d(PldO B )ij d{OldO B )ij 

l2 /v-.^ ,n vv! vvi f ^ w \ A t /~~\ A2 



(AA)jj —(Ab)u 
-{Aa)u (A B )jj 

where we used that (A A )jj(AB)jj = Aj (no summations), and the A^ are the eigenvalues 
of the matrix C = AB T . The remaining vN matrix elements below the block diagonal of 
dB T give a diagonal sub-matrix with v elements (A B )jj each. The resulting contribution 
to the Jacobian is 

/ N 

j= n i a '- a ?i ni^i"- ( 3 - 4 °) 

l<i<j<N j=l 

The prime on the product symbol denotes that only those factors with indices 
strictly below the block diagonal are to be included. Finally, we observe that the second 
product can be written as 

n ( n^ildet^r for even A^, 

IlK A s);;r = i ' ' (3.4i) 

U ( j li 1)/2 \ detB j\ U \( A B)NN\ u for odd N, 
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in which Bj is the j-ih 2x2 block along the diagonal of the matrix A#. The same 
statement is true in the more general case when (and A#) are not diagonal, as 



shown in Appendix A 



Writing everything together, we have for the total measure 
dAdB e v-TrAB^-±Tr(AAT + BBT) 



l N 



~ dA A dA B d/A A d/A B A l dO A B l dO B \\ l A ?- A ?in e?MV? ( 3 ' 42 ) 

l<i<j<N i=l 

e -^Tr(A5A A+ A5A B+ A5A A+ A5A fl ) { llf=i I det Bj\ v for even N 

\ Ili=r 1)/2 I det^HAfil^v for odd N. 

All constant factors are omitted here; we give the overall normalisation constant later. 
We reiterate that Bj here is the j-th 2x2 block on the diagonal of A B . The integration 
over the orthogonal dof as well as over the upper block triangular matrices Aa,b can 
now be performed as they decouple. The relevant dof for the right-hand side of eq. 
( I3.42p can thus be written in terms of the A^ and 2x2 blocks of the matrices A AjB 

( U^{^dB i e^-^+^e-^^ + ^\detB i \ v } for even N, 

i<] { n!fr 1)/2 { ditto }• dadb e v-A% e - v -±(aHi?)\ b \» for odd N. 

(3.43) 

We have thus reduced the problem of computing the jpdf to a simpler problem involving 
only 2x2 and lxl blocks, which can be handled just as the N = 2 and N — 1 cases 
that we treated in the previous two subsections. For this, it is necessary to order the 
eigenvalues as described in eq. (I2.10j) . Then, following eq. (13.321) . each 2x2 block will 
make the following contribution in variables Af = Zi (with A% i _ l and A^ either both 
real or complex conjugates, and ordered as described at the end of Section I3T21) 

d 2 z 2i -i d 2 z 2i (z 2i -i - z 2 i)w{z 2 i-i)w{z 2i ) 

x (%2i-i)%2i)©(z 2 i-i - x 2i ) - 2i5 2 (z 2i _ 1 - 4)e0/2i_o) . (3.44) 

In a similar way, we have for the lxl block when N is odd 

d 2 z N \z N \ v ' 2 e« g(z N ) 5{y N ) = d 2 z N w(z N ) 5(y N ) , (3.45) 

where z N will be real. Note that collecting all these quantities in eq. (13.431) for each 
i — 1 to [N/2] the factors A^_ x — A^ will combine with the Yl i<: AA. 2 — A 2 ) to make a true 
Vandermonde determinant ni<j( A f ~ A j)' an d we were able to drop the modulus sign 
because of the chosen ordering. The final answer for the jpdf is therefore as claimed in 
eq. ( 12.91) . The normalisation constant can be determined by keeping track of all volume 
factors and /^-dependencies; it is given by 

c N = (VO(N)2- N (2n)-^yy(2n)-^ 2 ^ 

(3.46) 
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4. Finite- and large-iV density correlation functions 

4-1- The kernel 

The kernel Kn{z\, z 2 ) as it initially appears in eq. (12 . 1 8[) is defined as follows [7j 

N N 

JC N ( Zl ,z 2 ) e^^^ 1 ^ 1 (4.1) 

k=l 1=1 

where the matrix A of dressed moments is related to T by 

A H = J d 2 Zl d 2 z 2 T{z x , z 2 ) z\- x z l 2 \ (4.2) 

However, to evaluate the matrix A and its inverse directly is not trivial in general. 
Fortunately, the kernel may also be derived from the expectation value of the product 
of two characteristic polynomials [20] 

H N (X, 7 ) = (det(A - V) det( 7 - V T )) N , (4.3) 

where T> was given in eq. ( 12.2ft . An explicit form for if at (A, 7) was derived by the authors 
in [20] Q In the case of the real Ginibre ensemble the spectral density was known prior 
to its integrability [4] and thus was used for the determination of the kernel in [7]. 

To establish the relationship between ]Cn(zi, z 2 ) and Hn_ 2 (\, 7), we first relate 
the latter to the complex eigenvalue density Rf N (z). If we choose 7 = A* with 
3mA 2 > 3m7 2 we have 

JV-2 



H N - 2 (\,\*)= [ d 2 Zl ... [ rf 2 ^_ 2 P J v_ 2 (^ 1 ,...,^_ 2 )(AA*)^r[(A 2 -^)(A 



1 -4) 



CN-2 1 



c N exp[-r/_(A 2 + A 2 *)]^(A 2 ) 2 (A 2 - A 2 *)(-2i)9(9mA 2 ) 

x / d 2 Zl ... d 2 z N _ 2 P N (z 1 , z N _ 2 , A 2 , A 2 *) , (4.4) 
Jc Jc 

where Pjv indicates that this jpdf is conditioned that the last two eigenvalues are complex 
conjugates and ordered. The last line in eq. (I4.4p is thus nothing but the complex density 
Rf N (\ 2 ), obtained by inserting a delta function into the partition function together with 
the constraint that the last two eigenvalues are complex conjugates. We thus arrive at 
the following relationship 

Rl N (z) = J W*M*') { _ 2i){z _ ,*) h n _ 2{V - Z] VF). (4.5) 

C N -2 \A 

On the other hand we have from eq. (12.251) an equation that relates the complex density 
directly to the kernel 

Ri,n(z) — / d 2 ulC N (z,u) T{u,z) (4.6) 
Jc 

+ Note that P and Q were called A and B in [20]; we set n in [20 to unity here. 
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and so for complex z inserting eq. (12.191) 

Rl N (z) = JC N (z,z*)w(z)w(z*)(-2i) sgn(y). (4.7) 
We can therefore make the identification (after analytic continuation in each argument) 

K N (u,v) = (u-v) , . 4.8 

cat_ 2 {uvyi 1 

Using the solution for H N (X,j) obtained in [20] as well as eq. (I3.46p . we then have for 
the properly normalised kernel that 

£n(u,v)= — — -— V — 77TTXT i Lj+i ( ~r~^~~ ) 



u 



8 7 r(4/xV) i/+1 \V+J ti + vV- I j+1 \Vw j \4(i 2 v- 

which is the same as eq. ( 12.201) . It can be further simplified to be expressed as an 
anti-symmetric derivative as follows. For modified Laguerre polynomials, we can use a 
recurrence relation to show that 



[j + 1) {L v j+1 (y)mx) - (x «- y)} = {y^L){y) + (j + v + 1 - y)L?(y)) L){x) 



-(x<->y) 
d_ _ 

dy dx 



y^--x^--(y-x))LUx)LUy). ( J 10) 



The kernel therefore becomes 



(4.11) 

where x and y are evaluated at x — zi/4fi 2 ri_ and y = z^j^i^f]- after the differentiations. 
The symmetric kernel in terms of Laguerre polynomials on the right hand side of eq. 
( 14. lip is nothing else but the kernel of the complex {(3 = 2) two-matrix model [16J. 

We have explicitly checked that a similar relation holds relating the kernel of the 
non-chiral real Ginibre ensemble (see the equation after 5.26 in second ref. [S]) to 
(d y — d x — 2(y — x)) operating on the kernel of the non-chiral complex Ginibre ensemble 
(see eq. 40 in [32]) and thus this is a more general feature. 

Note that the convention used in [2U] for defining the 'kernel' was different from 
that adopted here (in eq. (14.91) ); in that paper, there was no division by (uv) u ^ 2 , and 
the summation ran to N. There was also a minor typographical error in the arguments 
of one function (eq. 33 in [20]). 

4-2. Finite-N results 
From eq. (12.181) we have 

R x { Zl ) = f d 2 z)C N (z 1 ,z)^(z,z 1 ) = Rf(z l )+8(y 1 )Rf(x 1 ) . (4.12) 
Jc 
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We now simply insert the finite- ./V kernel eq. (14. 9 p and weight function eq. (12.121) into 
this using eq. (I2.19p . to give 

^ e -r,%t(z2+ z *Z)-± K „ (2 v lt\z\ 2 ) erfc (2 V+ Vi\%m, 



i a 

K 2l 

3=0 



V 



+ 1)! 



c.c. 



(4.13) 



V- 



87r(A/j, 2 r] 

N-2 
3=0 



rj_ 

V+ 



(/+| , dx'sgn(x - x , )\xx , \^ 2 e^ x ^2K^(ri + \x\)2K^(ri + \x , \) 

23 a + 1)! i . ,, / \ , „ / ./• 



L 



4/i 2 ?7_ 



J 1 ' 

L 3 



[x 



x) 



(4.14) 



These results are valid for even N only. Alternatively we could have used the form eq. 
( 14. lip which is more reminiscent of the corresponding chGOE result [35] at fi = 0. In 
fact we have checked that in the limit fi — > the complex density eq. (I4.13P vanishes, 
and the real density eq. (I4.14p reduces to the finite-iV expression (equation 5.18 in 
|35j ) for the chGOE, after using some identities for modified Laguerre polynomials. We 
have also checked these results numerically using Monte Carlo, by generating random 
matrices and explicitly diagonalising them. 

As the last step we can change from squared variables z = x + iy = A 2 to Dirac 
eigenvalues A, using eq. ( 12.28}) . These two densities are illustrated^ in Figures [T] and [2} 
showing the localisation of the support for finite N. 



-^1 Dirac( Z ) 




Figure 1. The complex spectral density i?f £>irac( z ) f° r finite N = 10 at maximal 
non-Hermiticity /i 2 = 1 (left) and intermediate /i 2 = 0.5 (right), both for v = 0. 
We show only the first quadrant for symmetry reasons. For /i = 1 we see a circular 
"support" growing with y/N , apart from the repulsion from the axes. For decreasing 
H the "support" becomes an ellipse, with the eigenvalues moving towards, as well as 
onto, the real axis. Note the increased height in the right plot. 



Here and in the following numerical integrals arc carried out using 36J. 
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Figure 2. The spectral density R^Dirac( x ) f° r rea ^ eigenvalues (blue full line) and 
R^Dirac(^ x ) f° r purely imaginary eigenvalues (red dashed line) for finite N = 10 at 
almost maximal non-Hermiticity /i 2 = 0.95 (left) and intermediate /i 2 = 0.5 (right), 
both for v = 0. Here we have chosen (almost) the same parameter values as Fig. [TJ 
with /i 2 = 0.95 here close to 1 there. Because of chiral symmetry real and imaginary 
eigenvalues come in ± pairs and we only have to show the positive axes, comparing 
both distributions in the same plot. Whilst for fi 2 — 0.95 the distributions of real 
and imaginary Dirac eigenvalues are almost the same, for /i 2 = 0.5 there are more 
eigenvalues on the real than on the imaginary axis. 



4-3. The large-N limit at strong non-Hermiticity 

In the strong non-Hermitian limit, we keep \i fixed as we take the large- N limit. This 
necessitates no rescaling of the eigenvalues (see Section EJ), which is why our result is 
also true away from the origin. Let us first determine the large- N limit of the kernel. 

Reference [31] eq. 8.976.1 gives us the so-called Hille-Hardy formula (the equivalent 
of Mehler's formula for Hermite polynomials): 

Si™*) - E - &^T*™I, (Mf ) . (4.15) 

We can insert this into eq. (14.111) and evaluate it after extending the sum to infinity On 
differentiating the right-hand side, certain terms cancel, and we can therefore establish 
that 

ay ox 1 — z 

Hence the limit as N — > oo of the kernel is easily seen to be 

JC^z 1 ,z 2 )= lim IC N (z 1 ,z 2 ) = ^(z 1 -z 2 )e^-^Hz 1 z 2 r^ 2 I,(2 V+ ^zTz- 2 ). (4.17) 

Because of this simplification this is now proportional to the (3 = 2 kernel at strong 
non-Hermiticity. Multiplication by the weight function eq. (12.121) which contains the 
modulus \ziz 2 \ u l 2 will only cancel the pre-factor in eq. (I4.17P up to a phase. Putting all 
ingredients together we can determine the eigenvalue densities using eq. (12. 25ft : 
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p% s (z) = sgn(Smz){-2i)(z - z*)'^I u {2 V+ \z\) 



, alt 
x 2 / — exp 



8tt 

2 / 2 i *2\± 1 



(4.18) 



Kil (2ri 2 + \z\H) erfc (Vf v 7 * | 



P^) = ^2i<: | (77 + |a;|) ^ dx'\x-x'\ 2K*( V+ \x'\) I„(2 V+ Vxl^) (4.19) 

+ f dx' \x - x'| 2K^(r] + \x'\) J v (2r] + ^f x\x'\) \ . 

J — oo / 

Note the change from Bessel-7 to Bessel-J function inside the integral for negative 
arguments, after taking into account the aforementioned phase. 

If we rescale the eigenvalues as 2r] + z 2 — > z 2 (and divide the densities by 2r] + 
accordingly) we obtain the same densities as at maximal non-Hermiticity. These are 
obtained by setting fj, — 1 =^ r] + — | above, or by starting from the kernel eq. (12.301) 
at maximal non-Hermiticity. This feature that the strong limit can be obtained by 
rescaling the case of maximal non-Hermiticity is generically true for complex RMT, see 

EDI. 



Dirac (A) 



P2 Dirac (A) 




Figure 3. The complex spectral density £j» rac (A = x + iy) at maximal non- 
Hermiticity fi = 1. It is shown only in the first quadrant for symmetry reasons, for 
v = (left) and v = 2 (right). Increasing the number of exact zero eigenvalues v leads 
to a stronger repulsion from the origin. At v — this repulsion is still present due to 
chiral symmetry (or technically speaking the presence of the Bessel-.fr function). 



In Figures [3] and H] we show the densities of Dirac eigenvalues A for complex 
eigenvalues and real or purely imaginary eigenvalues respectively, using the mapping 
eq. (I2.28p . Because of the rescaling property just mentioned we only show results 
here for maximal non-Hermiticity. One can check analytically using a saddle-point 
approximation including the fluctuations that for asymptotically large x, y the densities 
eqs. (14.181) and (14.191) decay as ~ l/|z| and l/y/x respectively. After the mapping to 
Dirac eigenvalues eq. (I2.28f) they thus reach a plateau as seen in the figures. We note 
that the profiles of the densities on the real and imaginary axis are very reminiscent to 
parallel cuts through the complex densities, for both values of v shown. 
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Figure 4. The real spectral density of Dirac eigenvalues on the positive half-line at 
maximal non-Hermiticity fi = 1 for v = (left) and v = 2 (right). Because of chiral 
symmetry it is symmetric on the negative real line, and because /i — 1 it is identical on 
the imaginary axis. For v = 2 we see the increased repulsion from the origin compared 
to v = 0, as for the complex eigenvalues in Fig. [3l 



4-4- The large-N limit at weak non-Hermiticity 



In the weak case, we scale /i and z with N as follows (as compared with the unsealed 
case, i.e. we implicitly include here the first rescaling discussed in Section 12.41) 

a 



(JL = 



<2N 



ANz 



(4.20) 



where a and z are kept fixed throughout. Because of the rescaling of the eigenvalues 
here we are magnifying the region around the origin. 
For the weight function we thus simply obtain 

\z\ 



lim K« 77+ 



2 



and 



lim exp 



Z ' 




" z ' 




= exp 


Ua 2 \ 



(4.21) 



For the kernel we are interested in the limit of eq. (14. 9p in terms of the rescaled variables, 
i.e. /j, r/_ and r/ + which all now depend on N. Instead of using eq. (14.111) it is slightly 
simpler if we can rewrite eq. ( 14. 9 p so that the Laguerre polynomials inside the sum are 
of the same degree j. Using the recurrence relationship for the Laguerre polynomials 
(eq. 8.971.2 in |34J ) , we have 

(n + 1) {U n+l {v) U n [u) -(„<->„)}=((„ + ! + v )U n {v) - vU n + \v)) U n [u) 



The kernel can therefore be written as 

JV-2 

" V- 



/Ctv [Zl, z 2) 



i V - 



2.7 



U + v) 



U *-* V ) 

(4.22) 
(4.23) 



T u 



Z-2 



4fi 2 r]- 



z 2 L» 



v+l 



Z2 



4fi 2 rj- 



L 3 



Zl 



4/i 2 r/_ 



We now wish to take the limit iV — > oo. For this to exist we must multiply the kernel by 
the spacing as well as by the appropriate number of zero-eigenvalues from the weight, 
as given in eq. (I4.27P below. In eq. (I4.23P we will replace the sum with an integral over 
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the variable t = -jfe E [0, 1]. Because of the different scaling in the weak limit we cannot 
use the Hille-Hardy formula as before. 

In detail, using eq. 8.978.2 in [34] we have for some real constant v and fixed t G [0, 1] 
the standard Bessel asymptotic of the modified Laguerre polynomials: 

'' X1 t v/2 x- v/2 J„(2v^t). (4.24) 



lim 



We also have (keeping t fixed so that j 

„2 \ 2i 



tN -> oo) 
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iV->oo \ 1 4- °il 
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r v . (4.25) 



'2N 



{\/^Jv+l{s^fzi)Jv{s^fz?) - (Zl Z 2 )} 



(4.26) 
(4.27) 
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z 2 -^- - zi-^- ) j als s e 2a s J u (sy^i) Ju {sy/z^ 



1287ra 2 

after a simple change of variables in the integral. In the last line we expressed the kernel 
as a derivative of the (3 = 2 kernel at weak non-Hermiticity. This corresponds to eq. 
(14. lip in which the second term — (y — x) now becomes sub-leading compared with the 
derivatives. The same relation between the large- N kernel at weak non-Hermiticity is 
true for the (3 = 1 and (3 = 2 Ginibre ensembles [H [32] respectively. 
Collecting all elements we have for the complex density 

1 



— 2i sgn(3m z) exp 



r 2 9ftez 



K w {z, z* 



(4.28) 
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In Figure [5] it is shown after mapping to Dirac eigenvalues. As a consistency check 
we can take the limit a — > oo while keeping z/a 2 fixed to obtain once more the complex 
density in the strong non-Hermiticity limit, eq. (|4.18|) . The precise mapping of weak to 
strong eigenvalues is given by 7% — > 2r/ + z. Whilst the matching of the integrals over t 
is straightforward the mapping of the kernels multiplied by the weight is more involved. 
Changing variables we have the following identity: 

• i 
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1 
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(4.29) 
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Figure 5. The complex spectral density £>$ roc (A = x + iy) at weak non-Hermiticity 
for parameters a 2 — 0.2 (left) and a 2 = 1 (right), both at v = plotted on the same 
scale. For increasing a more complex eigenvalues move in towards the imaginary axis, 
reaching Fig. [3] left in the limit a — ► oo. 

where in the last step we have extended the integral to infinity and used eq. 6.633.2 [34]. 
The last differentiation is trivial, in effect acting only on the exponential function. We 
thus obtain for the limiting kernel 

{W) • <"°> 

which precisely cancels the exponentials from the weight eq. (14.211) to arrive at the 
complex density at strong non-Hermiticity eq. (I4.18p . 

We now turn to the real density at weak non-Hermiticity. Looking at the definitions 
eqs. (14.131) and ( 14. 14ft the main difference to the complex density is that here the kernel 
is integrated, whereas the complex density is simply given by the kernel multiplied by 
the weight. 

Unfortunately, and in contrast to the strong non-Hermitian case, at weak non- 
Hermiticity the large- N limit and the integration do not commute, with the integral 
over the weak kernel (I4.27P not being absolutely convergent. Such a feature might have 
been expected, as the same phenomenon occurs for the chGOE at /i = [35]. However, 
in that case, the integrals could be computed exactly before taking the large- N limit, 
leading to the correct result, which differs from the naive limit by a factor of 2 in the 
normalisation^!- The integrals in eq. ( 14. 14ft at finite N are more involved, and this 
matter will be addressed in future work. 

For that reason we show in Fig. [6] the real density for finite but large N = 10 and 
20, using the weak scaling from eq. (14. 26ft and eq. (14. 20ft for a given a. This underlines 
that the weak limit for the real density does exist and convergence is rapid. Furthermore 
we have checked this by superimposing data for N = 100 using numerically generated 
random matrices. 

As a feature common to the strong limit we note that the densities of real and 
purely imaginary eigenvalues in Fig. [6] resemble cuts through the complex density in 
Fig. [5] left at the same value of a 2 = 0.2. 

jj A quick guess generalising this to our setting fails. 
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Figure 6. The real spectral density of Dirac eigenvalues on the positive half- line 
scaled as for weak non-Hermiticity for a 2 = /j, 2 /27V = 0.2 at v = 0. We show results 
for finite TV vs numerical simulations. Both the eigenvalue densities for real eigenvalues 
^? Dirad x ) an d f° r P ure imaginary ones £n rac {ix) are displayed in the same plot 
for comparison. Real density: N — 10 (blue full line), N = 20 (dark green dot-dashed 
line); imaginary density N — 10 (red dashed line) N — 20 (black, dotted line). Dots 
with error bars: N — 100 Monte Carlo simulation of 10 4 matrices. 



5. Conclusions 

In this paper we have solved the chiral extension of the Ginibre ensemble of real 
asymmetric matrices. It is given as a two-matrix model of rectangular matrices with 
real elements and depends on a non-Hermiticity parameter [i. This model is relevant 
for computing the non-Hermitian spectrum of Dirac operators with real elements in 
field theory. Our work completes the programme of solving the three chiral or Wishart- 
Laguerre counterparts of the classical Ginibre ensembles, where earlier works by Osborn 
and one of the authors extended the models with complex and quaternion real elements 
respectively. 

Whilst our model inherits most of the integrable structure of the real Ginibre 
ensemble its joint probability distribution required a more complicated calculation, 
which took much of our effort here. Just as in the Ginibre ensembles the probability 
density for the matrix elements is Gaussian, whereas the one for the eigenvalues becomes 
non-Gaussian. It contains a Bessel-.fr function and integral thereof, replacing the role 
of the complementary error function in the real Ginibre ensemble. 

The main building block for all eigenvalue correlation functions is given by a kernel 
of Laguerre polynomials in the complex plane and was derived in a previous paper. 
Here, we give all eigenvalue density correlation functions for finite (even) N valid for 
all values of /x, in particular the spectral one-point densities for real eigenvalues and 
for complex non-real eigenvalues. Moreover, we have uncovered a way of expressing the 
(3 = 1 kernel in terms of the (3 = 2 kernel at finite and large N, both for the chiral and 
non-chiral Ginibre ensembles. We conjecture that a similar relation holds for the (3 = 4 
kernel as well. 

When taking microscopic large-iV limits we focus on the origin where the chiral 
symmetry of our model is the most important. For both the limit at strong non- 
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Hermiticity with < /i 2 < 1, and the limit at weak non-Hermiticity with fi 2 ~ 

we give compact expressions for the kernel. This leads to explicit expressions for the 

complex spectral one-point densities and the real density at strong non-Hermiticity. 

It would be very interesting to compare these results with simulations from non- 
Hermitian lattice gauge theory, as was successfully done previously for the other two 
chiral two-matrix models. 

Further extensions would be to investigate the bulk or the soft edge scaling limit. 
We expect that in the former the chiral ensembles will agree with the Ginibre ensembles, 
as the effect of chirality becomes unimportant in the bulk. 

This work has been supported partly by European Network ENRAGE MRTN-CT- 
2004-005616 (G.A.), an EPSRC doctoral training grant (M.J.R) and the SFB/TR12 
of the Deutsche Forschungsgemeinschaft (H.-J.S.). We thank Tilo Wettig for useful 
discussions. 



Appendix A. Details on the calculation of the Jacobian 

In this appendix we give a few more details to complete the computation of the Jacobian 
from Section [3731 In particular we first give a precise ordering of matrix elements leading 
to a block-triangular Jacobian. Second we will treat the case of non-diagonal matrices 
A^ and A B . 

For the first purpose we repeat eq. f!3.37|) including matrix indices, after dropping 
the outside rotations: 

N N+u 

{dA)ij = J2(°AdO A )ik(^A + A A ) kj -J2( a a + A A ) ik {0 T B dO B ) kj 
k=i fc=i 

+ {dA A )ij + (dA A )ij , (A.l) 

N+u N 

(dB T ) lJ = J2(°l d °Bhk(A B + A B ) kj -^(A B + A B ) lk {0 T A dO A ) k] 

k=l k=l 

+ {dA B )ij + {dAsh ■ (A.2) 

We now give an ordering leading to a block-triangular Jacobi matrix, with variables 
{dA, dB T } in the columns and {dA A , dA B , dA A , dA B , A dO A , 0^dO B } in the rows. For 
the block diagonal matrices A AjB and upper block-diagonal matrices A AyB (of different 
size) this is trivial: we just group them together with the corresponding elements of dA 
and dB T . For example, when N is even and v > 0, this will give 

(dA) u , (dA) m (dA) 21 , (<L4) 22 , (dA) 33 , . . . , (dA) NN , (dB T ) u , . . . , (dB T ) NN , (dA) 13 , 
(dA) u , . . . , (dA) hN+I/ , (dA) 23} . . . , (dA) N , N+v , (dB T ) 13 , . . . , (dB T ) N _ 2 , N (A.3) 

versus 



(dA A ) n , (dA A ) 12 , (dA A ) 21 , (dA A ) 22 , (dA A ) 33 , . . , (dA A ) NN , (dA B ) n , . . , (dA B ) m , (dA A ) 13 , 
(dA A ) u , . . . , {dA A ) 1>N+u , (dA A ) 23 , (dA A ) N>N+v , (dA B ) 13 , . . . , (dA B ) N _ 2iN . (A.4) 
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The resulting sub-Jacobi matrix is clearly the identity matrix, and the order we have 
picked is arbitrary as long as we pair (dA)ij with (dA A )ij or (dA^)^-, and respectively 
for B. 

It remains to order the matrix elements dA and dB T below the block diagonal. In 
order to obtain sub-blocks as in eq. (13.391) we will always pair (dA)^ with the square 
part of (dB T )ij, with i > j and i, j G 1, . . . , N and finish with the rectangular part of 
(dB T )ij, with % = N + 1, . . . , N + v. For that we write the following partial differentials 
denoted by "|" from eqs. IjA.ljl and (TA~2l) 



p<j<i N+u 

(dA|) 4>i = (OldO A ) ip (A A ) PJ - ^A)i q (0 T B dO B ) q3 , (A.5) 
p=i q>i>j 

p<j<i N 

(dB T \) i>j = Y (0 T B dO B ) ip (A B ) pj - Y (&BUO T A dO A ) q3 . (A.6) 
p=i q>i>j 
Here we have already used that are upper triangular, and that i > j lets only the 
independent elements of the orthogonal differentials appear (we have chosen the below 
block diagonal ones). 

For the variables A AjB below the block diagonal, in order not to interfere with the 
block diagonals eq. (13.391) . we introduce the following ordering of elements: 

(hj) -< (i,p) if P < J , 

-< (q,j) if q> * ■ (A.7) 



Here -< implies that the matrix element 



has to appear before the one on the right qj; 



on the left, that is (dA)^ ((0^dO A )ij) 



Looking back to eqs. ( 1A.5j) and ( 1A.6I) 
this implies that the elements (dA)ij and (dB T )ij that depend on the most elements 
of (0^dO A )ij and (OgdO B )ij will appear first, leading to a lower triangular structure. 
However, the ordering eq. (1 A. 71) is not unique. We will proceed 2x2 block-wise (plus 
1x2 blocks for the last row when iV is odd) going down the diagonal, in order to 
preserve the block-diagonal structure when the matrices A AtB are not diagonal. 
We thus continue the labelling in eqs. (1A.3I) and (1A.4I) as 

(gL4) 32 , (dB T ) 32 , (dA) 31 , (dB T ) 3U (dA) 42 , (dB T ) A2} (dA) 41 , (dB T ) 4l , (dA) 54 , (dB T ) 54 , . . . , 
(dA) N1 , (dB T ) N1 (A.8) 

versus 

(OldO A ) w , (OldO B ) 32 , (OldO A ) 31 , (OldO B ) 31 , {0 T A dO A ) 42 , (OldO B ) 42 , {0 T A dO A ) 4l , 
(OldO B ) 41 , (0^dO A ) 54 , (OldO B ) 54} . . . , {0 T A dO A ) N1 , {0 T B dO B ) Nl . (A.9) 

It remains for us to order the rectangular part of (dB T )ij with i = N + 1, . . . , N + v 
vs the corresponding (0 B dO B )ij. Using the same order as in eq. (I A. 71) we complete our 
Jacobi matrix by 

(dB T ) N+1>N , (dB T ) N+ljN _ 1 , . . . , (dB T ) N+ i A , (dB T ) N+2i N, (dB T ) N+l , tl (A. 10) 

ff Our ordering is different from appendix A. 37 in [31j for the real Ginibre ensemble. 
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versus 



{0 T B dO B ) (0 T B dO B ) N+1,N-1, 

(0 T B dO B ) 



(0 T B dO B ) N+lil , (0 T B dO B , 



N+2,N, 



(A.11) 



In the second part of this appendix we will deal with the case of Aa.b being 2x2 
block matrices rather than diagonal, which was omitted in Section 13.31 We start with 
N even, and first only deal with the matrix elements dA and the square part of dB T 
below the block diagonal (eqs. ( 1A.8I) and ( 1A.9I) ). In this case it is no longer sufficient 
to study one pair of neighbouring elements as in eq. (I3.39p . but rather four pairs. This 
leads to the following 8x8 matrix with the order chosen above: In the columns we put 
(dA)ij +1 , (dB T ) itj+1 , (dA)ij, (dB T )ij, (dA) i+1 j + i, (dB T ) i+1 j +1 , (dA) i+1 j, (dB T ) i+ i tj , and 
in the rows the elements {JJ T A dO A)ij and (0 B dO B )ij in the corresponding order. Using 
eqs. ( 1A.1I) and ( 1A.2I) this leads to the following sub-matrices down the diagonal for each 
odd i and odd j: 
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(A.12) 

We can easily verify (using the symbolic manipulation capabilities of Mathematica |36j, 
for example) that the modulus of the determinant of this is identical to 



Ji 



'j 



(A 



SjDi) 



(A.13) 



in which Di is the determinant, and Si the trace, of the 2x2 matrix Ui given by (only 
relevant for odd i) 



(A.14) 



{\ A )u (A A )i,i+i \ f (A B )« (Ab)m+i 
(A A )i+i,i (A A ) 

The same definition obviously applies for subscript j. However, Di and Si can of course 
be written in terms of the eigenvalues of Ui, which we denoted A^ and Af +1 , i.e. 

Di = det Ui = klK 2 i+l , St = Tr U t = A 4 2 + k 2 l+l . (A. 15) 

Substituting these into eq. (1A.13j) . and factorising, gives 



Ji 



(A- -A?) (A 



2 

i+1 



A J 2 )(A l 2 -A| +1 )(A? +1 



(A.16) 



This is exactly what we got in the case when the A^ and A# matrices were diagonal, i.e. 
the product of four copies of eq. (13.391) . one for each of the combinations (i + 1, j), 

+ and + + 

When N is odd, we can treat the 2x2 blocks up to — 1 (inclusive) as before. We 
still have to consider the final column of dA, etc., and it is necessary here to treat (j, N) 
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and (j + 1, N) together as 2 x 1 blocks (for odd values of j in the range 1 < j < N). 
Hence, we have extra diagonal blocks in the Jacobian (for each odd j) of the form 

[ (^A)j+l,j+l — (Ab)atAT (^A)j+l,j \ 

— (^a)nN — (Ab)j + ij + i (AB)j+l,j 

\A~A)j,j+l (A A )jj —(AgjNN 

V (Aski+i -(Aa)nn (A B )jj J 

which is in fact the first sub-block of eq. (1A.12I) with i = N. The modulus of the 
determinant of this is identically equal to 



(A.17) 



Ji 



Nj 



D 3 - SjA% + A 



N 



(A.18) 



where we used that (Aa)nn(Ab)nn = Af^ (no sums), and Dj and Sj were defined as 
before. But we can switch to writing Dj and Sj in terms of the eigenvalues of Uj, and 
this gives 

,2 



Jl 



Nj 



(A^-A )(A 



4 



A 



■3+1' 



(A.19) 



We see that this is again of the expected form. 

Let us turn to the remaining variables, the rectangular part of dB T in eqs. ( 1A. 101) 
and (lA.lljl . We first suppose that is even. For each i and each odd j we find that 
the variables (dB T ) i)J+1 , (dB T )ij and (0^dOB)i,j+i, {O^dOs)^ are now coupled into the 
following 2x2 blocks: 



(A.20) 



(Ab) 3+1,3+1 (Ab) 3+1,3 
( A El)j,j+l (A B )jj 

whose determinant is (after a relabelling 2j — 1 — ► j) simply det Bj (using the definition 
of Bj after eq. (13.411) ). This is independent of i which takes v different values, 
Af + l<z<A^ + iA We therefore have a total contribution to the Jacobian of 

7V/2 

Y[ I det Bj\ u for N even, (A.21) 

3=1 

which is exactly as before (eq. (I3.4ip ). 

Finally, when A^ is odd, we just pick up an extra factor of (Ab)nn from each of the 
v cells in the last column (dB T ) i N 1 and so 

(7V-l)/2 



17 |det^n(AB) nn\ v for A^ odd. 



(A.22) 



This ends the calculation for non-diagonal matrices A^b- 
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